Density of states of a dissipative quantum dot coupled to a quantum wire 
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We examine the local density of states of an impurity level or a quantum dot coupled to a fractional 
quantum Hall edge, or to the end of a single one-dimensional Luttinger-liquid lead. Effects of an 
Ohmic dissipative bath are also taken into account. Using both analytical and numerical techniques 
we show that, in general, the density of states exhibits power-law frequency dependence near the 
Fermi energy. In a substantial region of the parameter space it simply reflects the behavior of the 
tunneling density of states at the end of a Luttinger-liquid, and is insensitive either to the value of 
the dot-lead interaction or to the strength of dissipation; otherwise it depends on these couplings too. 
This behavior should be contrasted with the thermodynamic properties of the level, in particular, 
its occupancy, which were previously shown to depend on the various interactions in the system only 
through the corresponding Fermi edge singularity exponent, and thus cannot display any Luttinger- 
liquid specific power-law. Hence, we can construct different models, some with and some without 
interactions in the wire (but with equal Fermi edge singularity exponents), which would have very 
different level densities of states, although they all result in the same level population vs. energy 
curves. 

PACS numbers: 73.23.Hk, 71.10.Pm, 73.20.Hb 



I. INTRODUCTION 

Understanding the behavior of low-dimensional elec- 
tronic systems has been one of the main challenges of 
experimental and theoretical physics in the last years. 
These systems are important not only as the basic build- 
ing blocks of nanoelectronic devices, but also for the in- 
tricate strongly-correlated phenomena they exhibit. An 
important subclass is that of metallic (gapless) one- 
dimcnsional systems, whose low energy dynamics is gov- 
erned not by Fermi liquid theory, but instead by the Lut- 
tinger liquid (LL) paradigmi. This description applies to 
a wide variety of experimental realizations, including nar- 
row quantum wires in semiconducting heterostructures, 
metallic nanowires, and carbon nanotubes. Closely re- 
lated are chiral LLs, formed at the edges of fractional 
quantum Hall effect (FQHE) systems^, and helical LLs, 
the edges of spin quantum Hall insulators'^. The effect of 
impurities on these systems is interesting from both the 
applicative and fundamental points of view. These impu- 
rities could also be intentionally introduced, in the form 
of, e.g., quantum dots and anti-dots. Hence, there is no 
wonder that such questions have attracted much effort 
recently. However, most of these studies were restricted 
to investigation of transport phenomenai^— , while other 
effects received much less attention^^r— . 

In this work we study probably the most basic exam- 
ple of such a system, namely, a single level in the vicinity 
of a fractional quantum Hall edge, or, equivalentlj»2S ^ a 
level attached to the end of a single LL wirei^. We will 
refer to the two components (in both systems) as "dot" 
and "lead" respectively. We include in our treatment 
the effects of short range dot-lead interaction, as well as 
the influence of an Ohmic dissipative bath (e.g., elec- 
tromagnetic fluctuations in gate electrodes In a 
recent work^^ we have studied the thermodynamic prop- 
erties of the model (e.g., the level population, entropy. 



and specific heat), and found that they are universal, 
in the sense that they depend on the various interac- 
tions in the model (intra-lead, dot-lead, and dot-bath) 
only through a single parameter, the Fermi edge singu- 
larity exponent. Thus, thermodynamics can neither be 
used to identify non- Fermi liquid behavior, nor to extract 
LL parameters. In this work we proceed to study, both 
analytically and numerically, the level density of states 
(LDoS), which may be probed by tunneling or absorp- 
tion spectroscopies. We find that the LDoS is sensitive 
to LL physics, even though its integral (times the Fermi 
function) gives the level occupancy, which is universal 
in the above sense. As we show below, the LDoS fea- 
tures power-law behavior near the Fermi energy. For not 
too strong interactions the exponent in this power-law 
is actually determined by LL physics alone, and is in- 
dependent of the level-lead and level-bath interactions. 
This and many other results derived below cannot be 
achieved using perturbative calculations^. 

The rest of this paper is organized as follows: In Sec.HIl 
we present our model, and apply to it the Anderson- Yuval 
Coulomb-gas (CG) expansion^S^— . We then proceed to 
analytic treatment of the LDoS in Sec. Illli and to nu- 
merical calculations in Sec. IIVI Finally, we summarize 
our findings in Sec. [V] 



II. MODEL AND COULOMB-GAS EXPANSION 



The system is described by the Hamiltonian H = 
Hd + Hl + Hdl + Hb + Hdb- The first term is the 
dot Hamiltonian Hd — eod}d, with d) and d the level 
creation and annihilation operators, respectively, and Eq 
the level energy. The second term is the lead Hamilto- 
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nian. It can be written in the form 
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(1) 



using chiral bosonic field (j){x) obeying the commutation 
relation [4>{x), (/>(y)] = i7rsgn(a; — y), where v is the veloc- 
ity of excitations^^. The level and the lead are coupled 
by: 

Hdl = tod^iO) + H.c. + Uo {dU - i) ^d^^iO). (2) 

The two terms in this equation describe, respectively, 
dot- lead hopping (with the tunneling matrix element), 
and local dot-lead interaction whose strength is Uq. The 
electronic annihilation operator at the end of the lead 
can be written as V'(O) = ;Ye*'^^°^/ v^/-\/27ra, where x is a 
Majorana operator, a is a short distance cutoff (e.g., the 
lattice spacing), and g is the LL interaction parameter 
{g < 1 for repulsion, g > I for attraction). For a FQHE 
system with filling g = v for electron tunneling (i.e., a 
dot outside the FQHE bar). Finally, the level is coupled 
to a bath of harmonic oscillator o^^'^^" — (describing, e.g., 
electromagnetic fluctuations in control gates), governed 



The dot-bath coupling can be 
\) Afc(aj,-|-afc). We assume 



by Hb = Y.k^ka\ak- 
written as Hdb = {d^ d - 

Ohmic dissipation, i.e., linear low- frequency behavior of 
the bath spectral function: Jb{^) = X^fe ^t^i^ ~ ^k) = 
Kuj. 

We examine this model employing the Anderson- Yuval 
CG expansion^!. In this approach, any quantity of inter- 
est is expanded to all orders in t^. This results in a series 
of correlation functions, which need to be evaluated for 
vanishing t^i^— The level- lead interaction gives rise to 
a potential at the end of the lead, which alternates be- 
tween Uq/2 and —Uq/2 whenever an electron tunnels in 
or out of the level. Similarly, the fcth bath oscillators ex- 
perience a shift in its equilibrium position, proportional 
to Afc . We thus have a sequence of Fermi edge singularity 
events^. The solution of this latter problem enables the 
calculation of all the terms in the series of correlation 
functions. 

RecentlySi we have studied in this way the partition 
function Z of the model, whose derivatives with respect 
to the parameters of the system (for example, the level 
energy eg and the temperature T) give us the thermo- 
dynamic properties (e.g., the level population, entropy, 
and specific heat). We were able to rewrite the series ex- 
pansion for Z in the form of a grand canonical partition 
function of a classical system of particles. These repre- 
sent hopping events generated by tg, and thus reside on 
the imaginary time axis of the original quantum model, 
which is a circle with circumference 1/T . Each particle is 
assigned a positive (negative) charge if it represents tun- 
neling of an electron from the lead to the dot (from the 
dot to the lead). Hence, there must be an even number 
of charges, which have to appear in alternating order of 
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FIG. 1: (Color online) A typical term (with 2A'' = 6 to-charges 
and s = — 1) in the CG expansions for (a) the partition func- 
tion [Eqs. ®-®]; (b) the dot Green function for t > t' 
[Eqs. ©-(IB])] (here s' = 1 and M = so a single to-charge 
precedes r', whereas there are 2M' = 4 fo-charges between r' 
and r). Signs and positions of the charges are indicated, with 
to-charges marked by thin green vertical lines, and d-charges 
by wide red ones. Level population as function of imaginary 
time is denoted by horizontal dashed blue lines. 



signs. The position of the ith particle is r^, and the sign 
of the charge of the first particle is denoted by s. The 
partition function then reads: 



Z = 



s=±l 
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rfT2 f dZl^^ScoisAn}) ^ (3) 



The charges have a fugacity y — yTo^/V, where Fq = 
ti'I^oPpl is the noninteracting level width [p^ — I/(7ru) 
is the corresponding lead local density of states], and 
^ ^ a/u is a short-time cutoff. The CG action is given 
by: 
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The first term of this classical Hamiltonian describes 
an interaction between the particles, with Vc{f) = 
ln{7rT^/ sin[7rT|f I]}. This interaction is similar in form 
to 2D Coulomb interaction, and is the origin of the name 
"CG expansion". The charges are two component vec- 
tors, where the two components correspond to the ef- 
fects of the coupling with the lead and the bath, re- 
spectively. They are given by Ci = s(— l)'~^eo, where 
the squared-magnitude of the charges, to be denoted by 
apES = |eoP, is the Fermi edge singularity exponent 
of the model. It is defined by behavior of the zero- 
temperature correlator of d^ip{0) with its Hcrmitian con- 
jugate, calculated at to = 0. This correlator decays as 
^-QFEs fQj. iQjjg time f. In our system we have found that 

6*0 = — 2y/g(5cff/7r, ^/K^ , where ^cff is the effective 

phase shift in the lead due to the dot- lead coupling^^. It 
is equal to Uq/{2v) in straightforward bosonization, but 
is more complicated in general. It may be extracted from, 
e.g., finite-size energy differences, which could be calcu- 
lated either numerically or analytically (via the Bethe 
ansatz)2i. The other part of the CG action accounts for 
the energetic cost of Eq per unit imaginary time for each 
interval in which the level is populated. Its form is analo- 
gous to an electric field applied to the classical system of 



charges. A typical configuration is depicted in Fig.[TJa). 

A similar treatment can be given to the LDoS, which 
we shall denote by pu^uj). It is equal to the imag- 
inary part of the level retarded Green function (mul- 
tiplied by — I/tt). The retarded Green function is in 
turn the result of analytic continuation of the Matsub- 
ara Green function from the upper half of the complex 
frequency plane'^'^. The latter Green function is defined 
by Gd(t - r') = ~TT{fre~"/^d{T)d'f{T')}/Z, where 
is the imaginary time ordering operator. Following the 
same methods as above, the numerator of this expres- 
sion can also be given a CG representation. This CG 
has the same form as Eqs. ©-(HI), with two additional 
charges of sizes ±ed, Cd = eb — (l/y^, 0), inserted at 
t' and r, respectively. These charges correspond to the 
level creation and annihilation operators appearing in the 
definition of the Green function. In the following we will 
refer to these as "d-charges" , to distinguish them form 
the other "io-charges" , which originate from the to term. 
The contribution of each such configuration is to be mul- 
tiplied by sgn(r' — r) to account for the Fermi statistics. 
Thus, for r > r' the full CG expression for the dot Green 
function is: 
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, (5) 



where s' = (1 — s)/2. The first 2M + s' tg-charges occupy the interval [0,r'], the following 2M' charges reside in the 
interval [t' ,t], and the last 2(iV — M — M') — s' tg-charges are in the interval [r, 1/T]. The classical action is given 
by: 

2N 2N 

ScgAs, r', {r,}) = |eop ^ s,SjVc{t, - n) -f gq • e<j ^ [Vc{n - r') - Vc{n - t)] - |e<j| Vc(t - t') + 

i<ij — l i—1 
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where the sign of the ith io-charge is Si = 
s(— l)*~^sgn(Tj — T')sgn(Ti — t). A typical configuration 
is shown in Fig.[TJb). Similar expressions hold for t < t'. 

Comparing the two CG expansions, the following ob- 
servation emerges: the CG expansion for the partition 
function contains only three parameters: Fq, So, and 
ctPES) while expansion for the Green function depends 
on g too (through Cd). Hence, the different interaction 



types (i.e., interactions in the wire, the dot- wire inter- 
action, and the dot-bath coupling) affect the partition 
function through a single parameter, the Fermi edge sin- 
gularity exponent apES- Thus, thermodynamic measure- 
ments cannot be used to distinguish between the different 
interaction types. In other words, one can construct very 
different models, whose interactions differ in strength and 
even in sign, which will have the same thermodynamic 
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properties, provided Fq, Eq, and apES are indeed the 
same'^^. On the other hand, the LDoS, which depends 
explicitly on g, wiU exhibit different behavior for these 
different systems. Hence, it can be used to extract the 
strength of intra-wire interactions, as we show below. 



III. ANALYSIS OF THE LEVEL DENSITY OF 
STATES 

As noted in our earher work^^, the CG obtained here is 
identical to the original Anderson- Yuval expansioi*21 for 
the anisotropic single-channel Kondo model^, demon- 
strating that the models are equivalent^^ii^. Under this 
mapping the level population becomes the magnetization 
of the Kondo spin (plus one half). Hence, Eq is analogous 
to a local magnetic field. Similarly, J^y is related to Fq, 
and Jz to ofpEs- The CG parameters obey the famous 
Kondo renormalization group (RG) equations^LM^ which 
read, in our notations: 



dy 
dlnC 

dlnC 



apES 



-Ay apES- 



(7) 

(8) 



Thus, the system considered can be in one of two phases, 
a strong coupling (antiferromagnetic-Kondo like) phase 
and a weak coupling (ferromagnetic-Kondo like) phase. 
The transition occurs, for small Fq, at apps = 2-1- 
0{y^Toa/v). In the weak-coupling phase the Coulomb 
charges form tightly-bound pairs. The level is thus 
effectively decoupled at low energies, resulting in its 
population being discontinuous as a function of eo at 
zero temperaturei^ii^ii^. In the strong-coupling phase 
free Coulomb charges proliferate. The impurity is well- 
coupled with the lead, so the level population is analytic 
in Eg J and could be extracted from the Bethe ansatz so- 
lution of the Kondo problem'^*. In particular, for small 
values of Sq one has n{eo) — 1/2 ^ —Eq/Tk, where 
Tk = (w/a)(Foa/i;)i/(2-"'='Es) jg the Kondo temperature 
(effective level width, reducing to Fq in the noninteract- 
ing case)i^. Hence, in this phase the population does not 
exhibit any nontrivial power-law dependence on Eq or T. 
The same applies to other thermodynamic quantities. 

What are the implications of this on the LDoS? As we 
now show, we typically find that at zero temperature we 



have a power-law behavior po (t^) 



in the vicinity of 



the Fermi energy, i.e., when |w| is much smaller than Tk 
in the strong-coupling phase, and than the bandwidth 
v/ain the weak-coupling phase. The values of 6 in the 
different regimes are summarized in Table HI It should 
be noted that when T > 0, or when the lead length L 
is finite, such power-law singularity will be smeared and 
become [max(|cj|, T, u/i)]''. We will now consider each 
phase separately. 



TABLE I: Summary of the analysis of Sec HITl In the yicinity 
of the Fermi energy (when \uj\ and T are much smaller than Tk 
in the strong-coupling phase, and than the bandwidth ^ v/a 
in the weak-coupling phase) we have poioj) ~ [max(r, li^l)]*, 
with the values of the exponent 5 in the different regimes 
denoted in the table below. When two values are given, the 
smaller one will give the dominant contribution^^. 



Phase 


1^1 < 


£o 


» £o 


Strong-coupling 




- 1 


1/ff-l 


Weak-coupling 


1/5- 


' 1 


a^ES - 1 or 1/9 " 1 



A. The strong-coupling phase 

Let us start from the strong-coupling phase. When |£o| 
is large enough (with respect to Tk)^ CG charges must 
appear in tightly-bound pairs, since large intra-pair sep- 
aration is suppressed by the level energy, as dictated by 
the last term of Eq. ©. The two d-charges added in 
the calculation of the level Green function will also be 
accompanied by two screening <o-charges for the same 
reason. The resulting configuration should thus resemble 
Fig.[2Ia). The leading contribution to the Green function 
will then come from the residual interaction of these par- 
tially screened d-charges, whose charges are (±1 times) 
eb — (Sd = (1/y^, 0). Thus, for large r — r' we have 
Gd{t - t') - sgn(r' - t)|t - t'\~^'3. This leads to 



the usual tunneling density of states 



singularity at the end of a LL wire^''*. 

When eo is small (with respect to Ik), renormalization 
effects are significant. In addition to the usual CG RG 
Eq. ([5]) representing the screening of interaction between 
to-charges by pairs of nearby io-charges (separated by 
^)— , one can write down similar equations for the flow of 
the Green function parameters^. It is easy to see that 
the coefficients of the logarithmic interaction between any 
two charges (either both io-charges, both d-charges, or a 
mixed pair) are renormalized in the same way by the 
pairs of nearby fo-charges, i.e.. 



d(e; 



rflnC 



= -2y2 (e^ • Co + eo • e^) , 



(9) 



where = 0,d. Thus, the combination |eo — = 
|eoP + |edp — 2eo ■ ed is invariant, and retains its initial 
value of 1/g. At the strong coupling fixed point y is 
large, so, by Eq. dH), |ej| — y^Oppg = 0, where asterisks 
denote fixed-point values. Thus, jepl = 1/y^, so that 

Gd{t' - r) - (r' - r)-l'^d*l' = (r' - t)-i/9, resulting 
again m pu (w) - Since this behavior holds at 

both large and small Eq values, it should also apply at all 
intermediate values. 

Further support for this result is obtained by analysis 
of the particular case apps = 1 (the Toulouse limiljSi), 
where the CG, and all thermodynamic properties, reduce 
to those of a noninteracting resonant level (for which 
g = 1, Uq = 0, and K = 0). A nontrivial (i.e., in- 
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teracting) realization of this condition, which still per- 
mits an exact calculation of the LDoS, is the case of 
no coupling to a bath {K — 0), but with g — 1/A 
and a corresponding compensating value of the dot- 
lead interaction. Then Cd — (—1,0) — — eo, so the d- 
charges have the same magnitudes as the corresponding 
to-charges, but the opposite signs. Comparison of the 
corresponding CG expansions thus shows that the level 
Green function of the interacting system is equal to a 
two-particle Green function of the noninteracting reso- 
nant level (fvV(0,T)d(T)d1'(r')V'^(0,r')), up to a factor 
of (27ra)~^sgn(T' — r). The operators at r and r' in this 
noninteracting two-particle Green function are similar to 
the to term in the Hamiltonian, but with 'ip{0) replaced 
by ip^{0) to account for the signs of the d-charges in the 
interacting system. After a straightforward evaluation of 
this two particle Green function by Wick's theorem, we 
find for the LDoS the following expression: 



n 



= 2/9Lacoth (^^^ Im 



1 _ £o - iTp 

2 2TTiT 



- 2eo + 2iro 
Lj — 2so 

- eg + iTg 
2mT 



(10) 



where ipiz) is the digamma function'*'' . For small \uj\ and 
T we indeed recover the [max(T, |aj|)]'* behavior appro- 
priate for g = 1/4, for all values of eg- 

Physically, the result is clear: in the strongly-coupled 
phase the level behaves, at low energies, as the last site 
of the lead, so its LDoS is similar to the local density of 
states near the end of a LL wirei*^, i.e., pd(w) ~ 
Interestingly, not only dot-lead interactions, but even 
coupling to the bath does not modify this behavior. As 
a result of this, the LDoS exhibits a power-law behavior 
with exponent which depends only on the LL parameter 
g and not on apES: i-e., on the interactions in the wire 
but not on the level-lead and level-bath couplings. This 
is in contrast with, e.g., the level occupancy, which de- 
pends on apES but not on g, as discussed above. Below 
we also test these predictions numerically. 



B. The weak-coupling phase 

We now turn our attention to the weak-coupling phase. 
Here all the tg-charges are bound in pairs, so the high |eo| 
results discussed above (which should also hold in this 
phase) actually carry over to low values of |eo|- It should 
however be remembered that then they compete with the 
contribution of the to = term in the CG expansion (the 
weak coupling fixed point), the term representing the in- 
teraction of the unscreened d-charges, which gives the 
LDoS a contribution of the form Ppair(w) ^ |w|"fbs-i^ 
with apgg = |edP- This is simply the LDoS of a tunnel- 
decoupled level, broadened from a delta peak to a power- 
law by the Anderson orthogonalities in the wire^ and in 
the Ohmic bath^iiS^. Note that since the level is effec- 
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FIG. 2: (Color online) Typical configuration of the CG ex- 
pansion for the dot Green function [Eqs. ©-([Gl); here r > r', 
s = -1, 2iV = 6, s' = 1, M = 0, and 2M' = 4, as in Fig.fHb)] 
for (a) large |eo| in both phases (and possibly small |eo| in the 
weak-coupling phased) — d-charges are screened; (b) small 
|eo| in the weak-coupling phase — d-charges are not screened. 
Notations are the same as in Fig. [1] See the text for further 
details. 



tively decoupled, the usual |a;|^/^~^ behavior at the end 
of a LL wire need not apply anymore. 

One could actually proceed to study higher order terms 
in the weak-coupling regime (and similarly, for large |eo| 
in the strong-coupling phase). For small tg the lead- 
ing correction is dressing of the above-mentioned charge 
configurations by a series of pairs of close-by io-charges 
(close- by since |eo| and/or apES are large), as depicted 
in Fig. [2] One then has to sum over all the terms sim- 
ilar to Fig. [2Ja) for large |eo| (i-e., for <C |£o| in the 
weak-coupling phase) or all the terms similar to Fig.[21[b) 
for small |eo| (i-e., for |a;| S> |eo|)- Since each pair 
has a very small dipole moment (due to the proximity 
of the charges), inter-pair interactions are negligible in 
a first approximation. This is actually an imaginary 
time variant of the noninteracting blip approximation 
(NIBA)^^'^-. The following argument can spare us the 
need of explicit calculations. For C/q = and K = 
the d-charges are noninteracting, = (0,0). The sum 
over all the terms with pairs of nearby tg charges would 
be the same (in the current approximation) as the Green 
function of a noninteracting system consisting of a level 
tunnel-coupled to a bath of noninteracting fermions with 



power- law local density of states pl(w) 



lapES — 1 



(with 



some appropriate high-energy cutoff and normalization) . 
For the latter system the Green function can be easily 
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evaluated to give G^jj{iuj) = [iuj — £o ~ : where 

the dot self energy is'^'^: 

oo 

s"^M^(to)' / ^^dn, (11) 



so that E5)(jw) ^ a;«FEs-i^ small |eo| we indeed see 
that YPjj{u]) is subdominant with respect to the nonin- 
teracting contribution only if apES > 2, which is exactly 
the condition for the weak-coupling phase for small io- 
Then, exactly at Eq = 0, a delta-function term appears 
at cli = in the expression for the LDoS (similarly to the 
situation at = 0), whose coefficient is determined by 
the requirement that the integral of the entire expression 
for the LDoS corresponding to G%{iuj) is unitjJ^. 

Before discussing nonzero Uq and K, it should be re- 
marked that this NIBA-like approximation exactly repro- 
duces the perturbative (in the tunneling to) approach em- 
ployed in Ref. [Tol, and would lead to similar predictions 
for the behavior of the level population. Both approx- 
imations are justified only in the weak-coupling phase 
(or when |eo| is large), but not in the strong-coupling 
phase, where to grows under RG flow and thus cannot be 
treated perturbatively, similarly to the exchange J^y in 
the equivalent Kondo problem!^. In the strong-coupling 
regime perturbative results predict correctly the depen- 
dence of the Kondo temperature on to (again, just like in 
the Kondo modeP*) and the qualitative behavior of the 
LDoS for 1/2 < .9 < 1 and eo 7^ (at [/q = and if = 0), 
but deviate from our previous conclusions in many other 
respects. For example, as we discuss below, at nonzero 
Uo and K our NIBA calculations indicate that the ex- 
ponent in the power-law behavior of the LDoS at low 
energy may depend on these interactions too, in contrast 
with the situation in the strong-coupling phase, where 
the corresponding exponent depends only on the LL pa- 
rameter g, as shown above. Moreover, even for vanishing 
dot-lead and dot-bath interactions (the case treated in 
Ref. 1£) the NIBA/perturbative expression given above 
does not agree with our previous analysis of the strong- 
coupling phase: for (a) 1/2 < g < 1 and Eq = or (b) 
g > 1 and any eo, NIBA would suggest that the LDoS 
i.e., with the opposite exponent to 



varies as 



1-1/9 



the one appearing in the density of states at the end of a 
LL wire. Perturbative results would thus imply that the 
LDoS may be enhanced for 17 < 1 or suppressed for g > 1, 
which is clearly at odds with both our previous results 
and the behavior of density of states at the lead edge. 
Moreover, integrating the perturbative LDoS leads to the 
prediction that the level population may have a power- 
law dependence on So at the strong-coupling phase (for 
1/2 < g < 2/3}JS, which is in contrast with the analytical 
behavior expected from the exact mapping of our model 
onto the Kondo problem, as discussed above. To summa- 
rize, NIBA/perturbative (in to) expressions do not hold 
in general in the strong-couplin g ph ase. It may be noted 
that the numerical data of Ref. [To! does not cover these 
regimes of the strong coupling phase for which perturba- 



tive calculations disagree with our previous analysis. 

Returning to the discussion of the NIBA approxima- 
tion for the weak-coupling phase, we will now treat the 
more general case, i.e.. nonzero Uo and K (which was 
not addressed in Ref. [10). Again, screened d-charges 
terms [Fig. UKa)] are dominant for \uj\ ^ |eo|, whereas 
unscreened d-charges terms [Fig. [2Kb)] are dominant for 
\uj\ |eo|'^ Let us start from the latter case. Now 
that there is interaction between d-charges, the CG ex- 
pression for the Green function Gd{t — t') contains an 
additional a factor of the form ^ \t ~ t'|^"fes (interac- 
tion of d-charges with the pairs of close-by ^o-charges is 
negligible). Turning to the frequency domain, the LDoS 
of the unscreened d-charges contribution [Fig. [IJb)] will 
thus be the convolution of the LDoS p%icd) associated 
with G'^j~i{iui) from the previous paragraph, with a func- 
tion Ppair(i^) ^ [wl^FEs-i^ which is simply the LDoS of a 
decoupled {to — 0) level as discussed above, i.e.. 



P£,(a;) = 2sgn(a;) / Pc(w - Jl)ppair(r2)df^ 



(12) 



at T = 0. Thus, the |a;|"FEs-i behavior of the decoupled 
level will survive for vanishing |eo|, due to the delta peak 
in G'jj{iuj). This behavior actually applies to the entire 
\uj\ ^ [eoI region, which is exactly where the contribution 
of terms similar to Fig. [DJb) is important. Higher order 
corrections will give extra powers of which are 

subleading since apES > 2. 

Similar considerations apply to the screened d-charges 
contribution [Fig. [UJa)] when |a;| ^ [eoj- For Uo ^ 
and/or K ^ 0, two corrections are due. The first correc- 
tion takes into account the factors coming from the in- 
teraction between each d-charge and the neighboring to- 
charge, which are power-laws in the time domain. Since 
the d-to charges form tightly-bound pairs, we can take 
the limits of integration over their separation to infinity. 
They thus yield factors of 



^leol^dr ^ r(l_ 



^ i\eo\0 



l — eo-Sd 



(13) 



where T{z) is the gamma function'^^, and eo-ej, — (c^fes + 
ttPES ~ l/5)/2- The correction is then the ratio between 
this expression and its value at = 0. Apart from this 
constant factor, one must compensate for the fact that 
the inter-pair interaction of these two d-^o pairs gives the 
Green function a factor which varies as |t — r'|~^/^ (since 
jed — ebj^ — 1/5)) instead of the |r — t'|^"'''=^ dependence 
used in the calculation of G^. Hence, p%{u!) should be 
convoluted here with |u;|i/ff-"FEs-i^ ppj. |^| ^ 
the LDoS will thus retain the behavior of the 

strong-coupling phase. 

To conclude the discussion of this NIBA-type approx- 
imation, a general expression for the LDoS which inter- 
polates between large and small [w/eqI limits^ is of the 



(a) (b) 




(c) (d) 




FIG. 3: (Color online) MC results for the LDoS as a function of frequency (measured with respect to the Fermi energy) for 
three different models, all with opes = 2/3 and r/Fo = 0.1. The three curves in each panel correspond to three different values 
of g: g = 1 (black pluses), g > 1 (green X's), and g < 1 (red circles), as indicated in the legends. The insets present the the 
LDoS at the Fermi energy as a function of the temperature by symbols (with the same code as before), together with a cyan 
line showing best-fit to the expected power-law behavior Pd{0) ~ y^/s^i (only the prefactor is fitted). The different panels 
correspond to different values of eo and K: (a) eo = and K — 0; (b) £o = ZFo and K — 0; (c) eo = and K = 1/6; (d) 
£o ~ SFo and K — 1/6. In each case the value of Uo was chosen according to the requirement apES = 2/3. 



form of Eq. (|T2|) . but with Ppair(^^) replaced by 



^-1 
n 



(14) 



All the low-energy results of this section are summarized 
in Table n 



IV. NUMERICAL CALCULATIONS 

In this section we present the results of numerical cal- 
culations, verifying the conclusions of our previous anal- 
ysis, i.e., that the LDoS at low energies features a power- 
law behavior with the power determined by LL physics 
only (in the strong coupling phase) , although, as we have 



shown beforeSi, its integral (the level population) is uni- 
versal, and cannot be used to extract LL parameters. 

To calculate the LDoS we used classical Monte-Carlo 
(MC) simulations on the CG expansion of dot Green 
function^^. The MC update procedure used is similar 
to the one employed recently for the closely-related con- 
tinuous time quantum MC algorithm^^. After obtain- 
ing the imaginary-time Green function it was Fourier- 
transformed to Matsubara frequencies, followed by ana- 
lytic continuation to real frequencies using the Fade ap- 
proximant technique"^*. This yields the retarded Green 
function, whose imaginary part is proportional to the 
LDoS'^'^. Below we present data in the non-perturbative 
strong-coupling region, which confirms the results of our 
previous analysis. Actually, in the weak-coupling phase 
(which is accessible analytically through the NlBA-like 
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FIG. 4: (Color online) DMRG results for the average level 
occupancy n_D as a function of its energy, eo, for a discrete 
model of the lead [Eqs. (I15|) - H16|) ]. The curves on which the 
symbols reside (which serve as guides to the eye) correspond 
to the various ofes values, where the larger apES the nar- 
rower the curve and vice versa. The widest curve has similar 
parameters to those used in the MC simulations presented in 
Fig. [3] On each curve there are three different choices of U 
and Ud (all giving the same ofes value), denoted by the three 
different symbol types. See the text for further details. 



approximation) MC simulations are not efficient, since 
there CG charges are rare and averaging very slow. In 
this sense, our analysis and numerical calculations are 
complementary. 

The results presented in the different panels Fig. [3] The 
values of Uq and K are varied, in a way which keeps 
apES constant at a value of 2/3. Hence, the occupancies 
as functions of eg are the same, as we also verify below 
(actually, the CG representation would predict exactly 
identical occupations). The LDoS curves are, however, 
markedly different: depending on whether (7 > 1, g < 1, 
or g = 1, they have a maximum, a minimum, or no special 
feature near the Fermi energy, respectively. In the inset 
we demonstrate that in all cases the LDoS at the Fermi 
energy exhibits power-law dependence on temperature, 
Pz?(0) ^ T^/9~i^ as found in the previous section (cf. 
Table U). 

For the sake of completeness, we will repeat here some 
of our previous data on the level occupancy^l. Since 
the MC simulation is based on the CG, to have an in- 
dependent check of the universality of the level popu- 
lation we employed the density matrix renormalization 
group (DMRG)22, algorithm, using block-sizes of up to 
256. DMRG is also better suited to ground state cal- 
culations, and thus complements the necessarily finite- 
temperature MC in this respect. The model used is a 
half-filled tight binding chain with nearest-neighbor in- 



teractions. It is described by the Hamiltonian: 

JV-l 

(15) 

where c| (c^) is the electronic creation (annihilation) op- 
erator at the ith site of the wire [i = 1 . . .N, with 
L = iVa), and t and U are, respectively, nearest- 
neighbor hopping amplitude and interaction strength. 
The low energy physics of this model is known to be 
governed by LL theory for not too large interactions 
(i.e., \U\ < 2t),^ with g TT/[2cos-\-U/2t)] and 
v/{2at) = vr^l - {U/2ty/[2cos-\U/2t)]. The dot is 
still governed by the same iJ^i = eod^d, and is coupled 
to the lead through: 

Hdl = tdc\d + H.c. + Ud {dU - i) (4ci - i) (16) 

td, Ud are related to the corresponding parameters of 
the continuum version Eq. (HI by io — id\/o,i and 
Uq — Uda, a being the lattice spacing. We have pre- 
viously shown that boundary conformal field theory ar- 
guments and the Bethe ansatz solution yield that here 
5eff = tan-i({/rf/V4t2-C/2).2i 

The level population is plotted in Fig. |4] as a function 
of Eq. Different curves correspond to different values of 
apESj as indicated in the legend. On each such curve 
there are three types of symbols, denoting DMRG data 
on three different models: (i) g — 1 (i.e., U — 0) but 
nonzero Ud; (h) 5 7^ 1 (nonzero U) but Ud = 0; (iii) both 
U and Ud are nonzero. All the models are without cou- 
pling to the bath {K = 0). The values of U and Ud in each 
model were chosen so as to give the same value of apES 
for each curve. For model (iii) we used U = ±0.5t, with 
sign opposite to that of model (ii). In all cases we chose 
td to get To = 10~H and used N = 100v/{at) sites. The 
results clearly show that the occupancy is indeed univer- 
sal, depending only on apESj and not on the strengths or 
signs of the interactions {U and Ud)- It should be noted 
that the widest curve has similar parameters to those 
used in the MC simulations (Fig. IS]). 



V. CONCLUSIONS 

To summarize, we have studied, both analytically and 
numerically, the LDoS of a level coupled to a LL and to 
an Ohmic bath over the entire parameter space. We have 
found that in general it exhibits a power-law dependence 
at low energies. In large parts of the phase space this is 
just the power-law behavior of the tunneling density of 
states at the end of a LL wire. Thus, a measurement of 
the LDoS there can be used to extract the value of the 
LL interaction parameter g. In other regions it is also 
affected by level-lead and level-bath interactions. In any 
case the LDoS is explicitly sensitive to the value of the 
LL parameter g, although the LDoS determines the level 
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population, which was found before to be universalsl, 
and thus not to feature any LL-specific power-law. 

Acknowledgments 

We would like to thank Y. Weiss for his invaluable help 
with the DMRG calculations, A. Schiller for many useful 



suggestions, and Y. Gefen for discussions. M.G. is sup- 
ported by the Adams Foundation of the Israel Academy 
of Sciences and Humanities. Financial aid from the Israel 
Science Foundation (Grant 569/07) is gratefully acknowl- 
edged. 



^ A.O. Gogolin, A.A. Nersesyan, and A.M. Tsvelik, 
Bosonization and Strongly Correlated Systems, (Cam- 
bridge University Press, Cambridge, 1998); T. Giamarchi, 
Quantum Physics in One Dimension, (Oxford University 
Press, Oxford, 2003). 

^ For a review see: A.M. Chang, Rev. Mod. Phys. 75, 1449 
(2003), and references cited therein. 

^ For a review see: M. Konig, H. Buhmann, L.W. 
Molenkamp, T. Hughes, C.-X. Liu, X.-L. Qi, and S.-C. 
Zhang, J. Phys. Soc. Jpn. 77, 031007 (2008), and refer- 
ences cited therein. 

^ C.L. Kane and M.P.A. Fisher, Phys. Rev. Lett. 68, 1220 
(1992); Phys. Rev. B 46, R7268 (1992); 46, 15233 (1992). 

^ A. Furusaki and N. Nagaosa, Phys. Rev. B 47, 3827 (1993). 

^ A. Furusaki, Phys. Rev. B 57, 7141 (1998). 

0. M. Auslaender, A. Yacoby, R. de Picciotto, K.W. Bald- 
win, L.N. Pfeiffer, and K.W. West, Phys. Rev. Lett. 84, 
1764 (2000). 

® H.W.Ch. Postma, T. Teepen, Z. Yao, M. Grifoni, and C. 

Dekker, Science 293, 76 (2001). 
^ A. Komnik and A.O. Gogolin, Phys. Rev. Lett. 90, 246403 

(2003); Phys. Rev. B 68, 235323 (2003). 

Yu.V. Nazarov and L.I. Glazman, Phys. Rev. Lett. 91, 

126804 (2003). 

" D.G. Polyakov and I.V. Gornyi, Phys. Rev. B 68, 035421 
(2003). 

1. V. Lerner, V.I. Yudson, and I.V. Yurkevich, Phys. Rev. 
Lett. 100, 256805 (2008). 

M. Goldstein and R. Berkovits, Phys. Rev. Lett. 104, 
106403 (2010). 

" A. Furusaki and K.A. Matveev, Phys. Rev. Lett. 88, 
226404 (2002). 

M. Sade, Y. Weiss, M. Goldstein, and R. Berkovits, Phys. 

Rev. B 71, 153301 (2005); 
1^ K. Le Hur and M.-R. Li, Phys. Rev. B 72, 073305 (2005). 
" Y. Weiss, M. Sade, M. Goldstein, and R. Berkovits, Phys. 

Stat. Sol. (b) 243, 399 (2006); Y. Weiss, M. Goldstein, 

and R. Berkovits, J. Phys.: Condens. Matter 19, 086215 

(2007); 

Y. Weiss, M. Goldstein, and R. Berkovits, Phys. Rev. B 

75, 064209 (2007); 76, 024204 (2007); 77, 205128 (2008). 
^® P. Wachter, V. Meden, and K. Schonhammer, Phys. Rev. 

B 76, 125316 (2007). 
^° G.A. Fiete, W. Bishara, C. Nayak, Phys. Rev. Lett. 101, 

176801 (2008); Phys. Rev. B 82, 035301 (2010). 

M. Goldstein, Y. Weiss, and R. Berkovits, Europhys. Lett. 



86, 67012 (2009); Proceedings of FQMT '08, Physica E 
42, 610 (2010). 

F. Elste, D.R. Reichman, and A.J. Millis, Phys. Rev. B 
81, 205413 (2010). 

A non-chiral LL with a boundary can mapped onto a chiral 
LL by unfolding the decoupled (Bogolubov-transformed) 
right- and left-movers [M. Fabrizio and A.O. Gogolin, 
Phys. Rev. B 51, 17827 (1995)]. 

A.J. Leggett, S. Chakravarty, A.T. Dorsey, M.P.A. Fisher, 
A. Garg, and W. Zwerger, Rev. Mod. Phys. 59, 1 (1987). 
U. Weiss, Quantum Dissipative Systems (World Scientific, 
Singapore, 1999). 

A. Kameney and Y. Gefen, Phys. Rev. B 54, 5428 (1996); 
arXiyxond-mat/9708109 

P.W. Anderson and G. Yuval, Phys. Rev. Lett. 23, 89 

(1969) ; G. Yuval and P.W. Anderson, Phys. Rev. B 1, 1522 

(1970) ; P.W. Anderson, G. Yuval, and D.R. Hamann, ibid. 
1, 4464 (1970). 

K.D. Schotte and U. Schotte, Phys. Rev. B 4, 2228 (1971). 
P.B. Wiegmann and A.M. Finkelstein, Zh. Eksp. Teor. Fiz. 
75, 204 (1978) [Sov. Phys. JETP 48, 102 (1978)]. 
^° Q. Si, and G. KotUar, Phys. Rev. B 48, 13881 (1993). 
M. Fabrizio, A.O. Gogohn, and P. Nozieres, Phys. Rev. B 
51, 16088 (1995). 

P. Nozieres and C.T. De Dominicis, Phys. Rev. 178, 1097 
(1969). 

G. D. Mahan, Many-Particle Physics (Kluwer, New York, 
2000). 

A.C. Hewson, The Kondo Problem to Heavy Fermions 
(Cambridge University Press, Cambridge, 1993). 
When 1/g < Op^s (i-^-, in the presence of strong dot- 
lead and/or dot-bath interactions), interaction between d- 
charges is large enough to make the screened d-charges 
contribution [Fig. [JJ^a)] dominant even for small |eo|- Using 
similar methods one can show that then the LDoS will vary 



|l/9-l 



just hke for \uj\ <^ \eo\ 



M. Abramowitz and LA. Stegun, Handbook of Mathemat- 
ical Functions (Dover, New York, 1964). 
^'^ P. Werner, A. Comanac, L. de' Medici, M. Troyer, and A.J. 
Millis, Phys. Rev. Lett. 97, 076405 (2006); P. Werner and 
A.J. Millis, Phys. Rev. B 74, 155107 (2006). 
H.J. Vidberg and J.W. Serene, J. Low Temp. Phys. 29, 
179 (1977). 

^8 U. SchoUwock, Rev. Mod. Phys. 77, 259 (2005); K.A. Hall- 
berg, Adv. Phys. 55, 477 (2006). 



